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U(l) gauge fields are decomposed into a monopole and photon part across the phase transition from the 
confinement to the Coulomb phase. We analyze the leading Lyapunov exponents of such gauge field configurations 
on the lattice which are initialized by quantum Monte Carlo simulations. We observe that the monopole field 
carries the same Lyapunov exponent as the original U(l) field. Evidence is found that monopole fields stay chaotic 
in the continuum whereas the photon fields are regular. 



1. Monopole and photon part of U(l) 

We begin with a 4c? U(l) gauge theory 
described by the EucHdean action S'{?7/} = 
/3 I]p(l - coseip), where Ui = t/a;,,^ = exp(i6'a;,^) 
and dp = Ox.p. + O^+j^^y - O^+p,^ - 9^,^ . We are 
interested in the relationship between monopoles 
and classical chaos across the phase transition at 
j3c ~ 1.01. Following Ref. [1], we have factorized 
our gauge configurations into monopole and pho- 
ton fields. The U(l) plaquette angles 9x,^u are 
decomposed into the "physical" electromagnetic 
flux through the plaquette 9x,iiu and a number 



of Dirac strings through the plaquette 



2iT m 



(1) 



where 9x.^v G (^tT; +'"'] and mx^fj,u 7^ is called a 
Dirac plaquette. 

2. Classical chaotic dynamics from quan- 
tum Monte Carlo initial states 

Chaotic dynamics in general is characterized 
by the spectrum of Lyapunov exponents. These 
exponents, if they are positive, reflect an expo- 
nential divergence of initially adjacent configu- 
rations. In case of symmetries inherent in the 
Hamiltonian of the system there are correspond- 
ing zero values of these exponents. Finally neg- 
ative exponents belong to irrelevant directions 
in the phase space: perturbation components in 
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these directions die out exponentially. Pure gauge 
fields on the lattice show a characteristic Lya- 
punov spectrum consisting of one third of each 
kind of exponents [2]. Assuming this general 
structure of the Lyapunov spectrum we inves- 
tigate presently its magnitude only, namely the 
maximal value of the Lyapunov exponent, imax- 
The general definition of the Lyapunov expo- 
nent is based on a distance measure d{t) in phase 
space. 
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In case of conservative dynamics the sum of all 
Lyapunov exponents is zero according to Liou- 
ville's theorem, ^ = 0. We utilize the gauge 
invariant distance measure consisting of the lo- 
cal differences of energy densities between two 3d 
field configurations on the lattice: 



d := — V ItrC/p - tiU'p 
Np ' ^ 



(3) 



Here the symbol stands for the sum over all 
Np plaquettes, so this distance is bound in the 
interval (0, 2iV) for the group SU(N). Up and Up 
are the familiar plaquette variables, constructed 
from the basic link variables Ux i, 



Ux,^ = exp UAl {T^ 



(4) 



located on lattice links pointing from the posi- 
tion X = {xi,X2,x^) to x -I- aci. The gener- 
ators of the group are — —igT'^/2 with r'^ 
being the Pauli matrices in case of SU(2) and 
j is the vector potential. The elementary 
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Figure 1. Exponentially diverging distance in real time of initially adjacent U(l) field configurations on 
a 12^ lattice prepared at /3 = 0.9 in the confinement (left) and at /3 = 1.1 in the Coulomb phase (right). 



plaquette variable is constructed for a plaquctte 
with a corner at x and lying in the ij-planc as 
Ux,ij = Ux,iUx+i,jUl_^_j ^Ulj. It is related to the 
magnetic field strength j.: 



U. 



exp (eijkaBl uT'' 



(5) 

The electric field strength E"^ • is related to the 
canonically conjugate momentum P^^i = Ux,i via 
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The Hamiltonian of the lattice gauge field sys- 
tem can be casted into the form 
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{U,V) 
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Here the scalar product stands for (A, B) = 
^tr(>liJt). The staple variable V is a sum of 
triple products of elementary link variables clos- 
ing a plaquctte with the chosen link U . This way 
the Hamiltonian is formally written as a sum over 
link contributions and V plays the role of the clas- 
sical force acting on the link variable U . 

We prepare the initial field configurations from 
a standard four dimensional Euclidean Monte 
Carlo program on a 12^ x 4 lattice varying the 
inverse gauge coupling /? [3] . We relate such four 
dimensional Euclidean lattice field configurations 
to Minkowskian momenta and fields for the three 
dimensional Hamiltonian simulation by identify- 
ing a fixed time slice of the four dimensional lat- 
tice. 



3. Chaos, confinement and continuum 

We start the presentation of our results with 
a characteristic example of the time evolution of 
the distance between initially adjacent configura- 
tions. An initial state prepared by a standard four 
dimensional Monte Carlo simulation is evolved 
according to the classical Hamiltonian dynamics 
in real time. Afterwards this initial state is ro- 
tated locally by group elements which are chosen 
randomly near to the unity. The time evolution of 
this slightly rotated configuration is then pursued 
and finally the distance between these two evolu- 
tions is calculated at the corresponding times. A 
typical exponential rise of this distance followed 
by a saturation can be inspected in Fig. 1 from 
an example of U(l) gauge theory in the confine- 
ment phase and in the Coulomb phase. While 
the saturation is an artifact of the compact dis- 
tance measure of the lattice, the exponential rise 
(the linear rise of the logarithm) can be used for 
the determination of the leading Lyapunov expo- 
nent. The left plot exhibits that in the confine- 
ment phase the original field and its monopole 
part have similar Lyapunov exponents whereas 
the photon part has a smaller Lmax- The right 
plot in the Coulomb phase suggests that all slopes 
and consequently the Lyapunov exponents of all 
fields decrease substantially. 

The main result of the present study is the de- 
pendence of the leading Lyapunov exponent Lmax 
on the inverse coupling strength /3, displayed in 
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Figure 2. Lyapunov exponent of the decomposed 
U(l) fields as a function of coupling. 

Fig. 2. As expected the strong coupling phase 
is more chaotic. The transition reflects the crit- 
ical coupling to the Coulomb phase. The plot 
shows that the monopole fields carry Lyapunov 
exponents of nearly the same size as the full U(l) 
fields. The photon fields yield a non-vanishing 
value in the confinement ascending toward (3 = 
for randomized fields which indicates that the de- 
composition works well only around the phase 
transition. 

An interesting result concerning the continuum 
limit can be viewed from Fig. 3 which shows the 
energy dependence of the Lyapunov exponents for 
the U(l) theory and its components. One ob- 
serves an approximately linear relation for the 
monopole part while a quadratic relation is sug- 
gested for the photon part in the weak coupling 
regime. From scaling arguments one expects a 
functional relationship between the Lyapunov ex- 
ponent and the energy [2,4] 



L{a) (xa''-^E''{a), 



(8) 



with the exponent k being crucial for the contin- 
uum limit of the classical field theory. A value 
of fc < 1 leads to a divergent Lyapunov expo- 
nent, while k > 1 yields a vanishing L in the 
continuum. The case fc = 1 is special leading to 
a finite non-zero Lyapunov exponent. Our analy- 
sis of the scaling relation (8) gives evidence, that 
the classical compact U(l) lattice gauge theory 
and especially the photon field have w 2 and 
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Figure 3. Comparison of average maximal Lya- 
punov exponents as a function of the scaled av- 
erage energy per plaquette ag^E. The full U(l) 
theory shows an approximately quadratic behav- 
ior in the weak coupling regime which is more 
pronounced for the photon field. The monopole 
field indicates a linear relation. 

with L{a) ^ a regular contimmm theory. The 
monopole field signals 1 and stays chaotic 
approaching the continuum. 
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